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Abstract  A  study  is  conducted  on  a  behavioral  model  to  study  man-machine  systems  when 
exposed  to  physiological  stress.  Using  such  a  model,  an  index  related  to  human  information 
processing  is  developed.  By  applying  this  model  to  empirical  data  from  a  centrifuge 
experiment,  the  index  is  calculated  ami  shows  some  interesting  results.  By  using  such  an 
index  in  conjunction  with  this  model,  the  experimenter  can  develop  a  “design  rale”  on  how 
to  choose  tracking  performance  tasks  to  further  study  different  types  of  human  tracking 
strategies. 


1.  INTRODUCTION 

The  behavior  model  concept  [1]  is  an  interesting  approach  in  relating  types  of  behavioral 
responses  to  a  parametric  representation  of  a  model.  If  the  model  can  classify  different  human 
responses  into  some  parametric  representation,  then  a  variety  of  studies  are  possible.  One 
could,  for  human  experiments,  vary  an  independent  variable  such  as  an  environmental 
stressor  and  observe  the  changes  in  the  empirical  data  resulting  from  the  man-machine 
interaction.  Once  an  accurate  cause-effect  relationship  can  be  developed  from  the  environ¬ 
mental  stressor  to  human  tracking  performance,  the  model  can  then  be  used  to  predict 
changes  in  human  behavior  for  a  wide  variety  of  different  levels  of  the  environmental  stressor. 
This  becomes  an  aid  in  the  design  of  experiments  to  test  protective  equipment  or  other  devices 
to  enhance  the  human’s  tolerance  to  *  stress  environment.  It  is  necessary  in  the  design  of  these 
experiments  to  not  only  design  the  stress  level  to  be  of  sufficient  magnitude  to  incur  a 
physiological  stress  on  the  subject  but,  in  addition,  it  is  necessary  to  design  the  tracking  task 
to  provide  a  mental  workload  requirement  that  is  sufficiently  challenging. 


O  2.  A  DESCRIPTION  OF  THE  MAN-MACHINE  INTERACTION 

♦  CJ 

Figure  1  illustrates  a  Mock  diagram  description  of  the  man-machine  experiment  that  will 
LtJ  be  considered  in  this  paper./(f)  is  a  variable  which  describes  the  target  forcing  function  which 
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Fig.  1.  Open-loop  stress. 


appears  on  a  cathode  ray  oscilloscope.  The  human  observes  on  the  display  the  closed  loop 
tracking  error  [dentoed  here  as  e(t)]  and  makes  stick  commands  which  become  the  input  into 
a  block  labeled  “machine.”  The  output  of  the  machine  is  the  variable  x(t)  which  is  compared 
to  the  reference  trajectory /(/)  resulting  in  the  error  e(t)  being  the  difference  between /(/)  and 
jc(/).  As  the  human  makes  stick  commands,  these  inputs  drive  the  cab  roll  dynamics  of  a  3 
degree-of-freedom  simulator  (centrifuge)  which  induces  an  acceleration  stress  denoted  as 
+  Gz.  The  term  Gz  means  an  acceleration  stress  vector  pointing  down  the  spine  of  a  subject. 
The  acceleration  stress  levels  created  by  the  centrifuge  produce  physiological  discomfort. 
Under  this  effect  the  subject  reduces  the  levels  of  his  stick  commands  and  does  not  attempt 
to  track  the  target  /(r)  with  the  accuracy  he  has  under  the  no  stress  condition.  In  th.s  sense 
the  subject  trades  off  “pain  to  performance.” 

3.  THE  BEHAVIORAL  MODEL 

An  important  aspect  of  human  tracking  involves  behavior  changes  [2].  Studies  of  human 
tracking  indicate  that  behavior  modifications  [which  ate  manifested  by  dramatic  changes  in 
characteristics  of  the  tracking  error  <>(/)]  can  also  be  due  to  properties  of  the  task  variable 
/(r).  This  is  best  illustrated  in  Fig.  2  in  which  a  dichotomy  of  behavior  can  now  be  explicitly 
stated. 

In  Fig.  2,  the  variable  /denotes  the  velocity  (time  rate  of  change  of  position)  of  the  task 
variable  f(t).  /  denotes  the  acceleration  of /(/)  with  respect  to  time.  From  previous  studies 
on  human  tracking  it  has  been  determined  that  for  small  values  of  /  and  /  (corresponding 
to  behavior  region  I  inside  the  circle),  the  tracking  task  is  easy  to  follow  and  the  subject 
maintains  high  levels  of  accuracy  in  tracking.  If,  however,  the  trajectory  f(l)  gets  outside  a 
critical  boundary  (and  enters  into  region  II — the  difficult  tracking  region),  a  behavior  change 
is  known  to  occur.  The  target  variable  /(/)  has  excessively  large  values  of  velocity  and 
acceleration  in  region  II,  which  results  in  the  human  exhibiting  “regressive  behavior.”  The 
term  regressive  behavior  is  borrowed  from  the  psychological  literature  which  means  that  the 
subject  will  substantially  reduce  his  stick  commands  and  not  try  to  track  the  target.  It  is  to 
his  benefit  (in  a  performance  sense)  to  exhibit  this  type  of  behavior  because  in  region  II, 
increasing  stick  commands  only  result  in  larger  error  signals. 

Thus  it  is  seen  that  human  tracking  behavior  defined  in  this  context  is  influenced  by  the 
following  two  variables: 

(1)  Physiological  stress— In  this  case  the  human  will  trade  off  tracking  performance  to 
discomfort  induced  by  the  environment. 

(2)  Characteristics  of  f(i) — If  the  task  moves  too  quickly,  regressive  behavior  is  known 
to  occur. 
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It  is  desired  to  incorporate  both  of  these  concepts  into  a  behavior  type  model.  To  achieve 
this  result.  Fig.  2  illustrates  this  concept.  The  behavior  model  assumes  the  following: 

(1)  A  critical  boundary  exists  which  separates  region  I  from  region  II.  For  simplicity,  this 
boundary  will  be  assumed  to  be  characterized  by  a  circle. 

(2)  Inside  this  boundary  defines  region  I,  an  easy  tracking  region  represented  by  small 
absolute  values  of  /  and  /  (inside  the  circle). 

(3)  Outside  this  boundary  defines  region  II,  where  the  human  behavior  is  characterized 
by  regressive  responses.  This  is  assumed  to  correspond  to  large  values  of  [/|  and/or  [/|. 

(4)  All  human  tracking  behavior  is  divided  into  only  two  regions.  This  assumption  seems 
tenable  based  on  discussions  with  pilots  involved  in  air  combat  situations  in  which  aircraft 
try  to  out  maneuver  each  other. 

(5)  The  critical  boundary  shrinks  or  expands  under  the  following  conditions: 

(a)  For  severe  levels  of  environmental  stress,  the  boundary  will  shrink.  This  indicates  a 
smaller  capable  region  of  tracking. 

(b)  The  addition  of  assistive  devices  (e.g.,  a  lead  angle  computer  for  a  gun  system)  may 
increase  the  size  of  the  boundary.  This  has  been  shown  to  be  the  case  in  certain  gun  systems 
[3]. 

(c)  A  good  tracker  will  have  a  larger  tracking  capability  region  as  compared  to  a  poor 
tracker.  This  has  also  been  shown  empirically  [2]. 

4.  A  STATISTICAL  ALGORITHM  TO  VALIDATE  THE  BOUNDARIES 

The  one  element  remaining  in  this  type  of  model  is  a  decision  rule  to  describe  in  some 
manner  the  critical  boundary.  In  Fig.  3,  a  plot  of  e(r)  versus  e(t)  is  constructed  using  similar 
arguments  as  described  previously  [2].  Region  I  can  be  assumed  to  be  within  an  elliptical  area; 
region  II  will  appear  outside  this  area.  It  is  also  shown  in  the  appendix  that  the  variable  e(t) 
is  actually  a  stochastic  process  e(t,  x)  but  is  denoted  here  in  the  text  of  this  paper  as  e(i). 


«( t, x )  r 
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The  distribution  functions  f(x)  and  /5(jt)  which  characterize  this  phase  plane  boundary  in 
Fig.  3  are  interrelated  by  a  correlation  relationship.  This  is  discussed  in  more  detail  in  the 
appendix. 

In  Fig.  4,  the  statistical  decision  rule  from  [2]  is  illustrated  in  a  heuristic  manner.  Data 
values  obtained  from  either  \e\  or  \e\  are  plotted  versus  time.  During  the  region  I  types  of 
behavior,  the  dependent  variable  (|ej  or  Jej)  has  a  value  close  to  a.  This  corresponds  to  the 
axis  of  the  variable  of  interest  in  the  ellipse  in  Fig.  3.  At  the  time  t,  (cf.  Figs.  3  and  4),  the 
value  of  this  variable  changes  from  approximately  at  op  units.  The  dependent  variable  must 
maintain  a  level  of  p  units  for  at  least  t2  sec  to  be  considered  into  region  II.  The  statistical 
decision  rule  depends  on  the  relative  values  of  p  and  a  and  is,  in  itself,  a  random  variable. 
Thus  there  may  be  some  error  in  defining  the  critical  boundary  (it  is  not  deterministic),  but 
this  error  can  be  obtained  from  the  properties  of  the  distribution  functions /,(*)  and /,(*). 

5.  LINKS  TO  INFORMATION  THEORY 

It  is  desired  in  this  paper  to  relate  the  model  thus  described  in  such  a  manner  analogous 
to  an  information  theoretic  method.  The  calculation  of  channel  capacity  in  information 
theory  is  no  simple  task.  Such  a  calculation  depends  on  properties  of  the  input  and  output 
signals,  the  noise  processes,  and  the  type  of  channel  considered  (continuous  or  discrete).  For 
the  special  case  of  a  bandlimited  input  signal  with  bandwidth  W  (rad/sec),  and  a  time 
continuous  Gaussian  channel  [4],  the  capacity  of  such  a  channel  can  be  written 

Cf^tV  log(l  +  S/N) ,  (1) 

where  S  and  N  are,  respectively,  the  signal  power  and  the  noise  power  in  the  channel.  The 
quantity  Cf  (units  of  bits/sec)  in  Eq.  (1)  provides  a  useful  bound  in  determining  actual 
transmission  rates  in  a  channel  and  also  in  calculating  capacities  for  other  types  of  channels 
(e  g-,  a  time-discrete  Gaussian  channel).  It  can  be  shown  [4J  that  C,  provides  an  upper  bound 
for  all  transmission  rates  in  a  channel. 

1*1  or  |*| 


Fig.  4.  The  fUdstkei  daemon  rate. 
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The  variable  of  interest  in  Eq.  (1)  is  the  ratio  (SIN).  In  the  field  of  communications,  the 
figure  of  merit  (S/N)  describes  the  important  aspect  of  a  channel  in  detail.  For  example,  if 
(S/N)  is  10  or  greater,  it  is  generally  accepted  as  a  good  channel  for  transmission  of 
information.  In  an  analogous  manner,  variables  similar  to  (S/N)  find  application  in  the  fields 
of  psychophysics  under  differerent  labels.  For  example,  in  the  area  of  manual  control, 
tracking  experiments  have  shown  for  hand  control  studies  what  is  referred  to  as  a  Weber’s 
law  effect.  If  Nx  represents  the  signal  power  in  the  hand  control  activity  not  correlated  to  a 
tracking  task  and  Sx  represents  the  total  signal  power,  the  ratio 

20  \ogm(NJSx)  =  —  20  dB  (2) 

is  a  well  known  relationship  determined  from  remnant  studies  [S,  6]  which  holds  for  the  low 
frequency  range  where  most  of  the  control  (stick  power)  is  known  to  occur.  The  relationship 
from  Eq.  (2)  implies  that 


SJNX  &  10  (3) 

holds  for  normal  neuromotor  tremor  response  in  this  type  of  situation.  A  similar  type  of  effect 
is  known  to  occur  in  situations  involving  the  lifting  of  a  weight.  If  a  human  were  to  lift  a 
weight  of  Wx  pounds  above  his  head,  and  if  A  Wx  was  the  root  mean  square  measure  of  the 
variability  (or  tremor)  associated  with  this  weight,  then  the  ratio 

(A  W,)/  =  constant  (4) 

is  known  to  hold  for  a  wide  range  of  H,  values  up  to  a  maximum.  This  is  sometimes  referred 
to  as  the  “Henneman’s  Size  Principle”  and  implies  that  variability  increases  in  proportion 
to  the  disturbing  force.  Thus,  for  larger  weights,  one  sees  larger  variability,  as  is  expected. 

A  fourth  area  where  this  type  of  information  theoretic  concept  occurs  is  in  parameter 
identification  procedures.  If  n  denotes  an  estimate  of  a  parameter  and  a  its  variance  in  this 
estimate,  the  ratio 


/  =  (5) 

becomes  a  valid  indication  of  system  order  [7J.  When  /  is  minimized  with  respect  to  the 
unknown  parameters  a  and  ji,  the  system’s  order  can  be  determined.  A  rigorous  test  on  system 
order  can  be  made  using  this  type  of  index.  If  /  exhibits  values  less  than  0.1,  the  order  of 
a  system  can  be  determined  within  accepted  scientific  levels  (p  <  0.0S).  An  interesting 
discussion  on  the  use  of  invariant  constants  can  be  found  in  the  classical  paper  [8]  by  Miller 
related  to  distinguishing  a  variety  of  psychophysical  stimulus.  Discussions  in  (9, 10]  also 
provide  motivation  for  studies  of  this  type  although  it  has  been  known  [1 1]  that  information 
theoretic  models  of  humans  involved  in  tracking  tasks  have  not  found  wide  application. 

The  ability  of  humans  to  act  like  information  channels  yields  some  provocative  results. 
If  capacity  is  measured  in  units  of  bits/sec,  studies  related  to  a  wide  number  of  reading 
experiments  yield  [14]  rates  on  the  order  of 40-50  bits/sec  for  this  type  of  information  channel. 
In  terms  of  vision,  however,  Kelly  [IS]  measured  the  information  capacity  of  a  single  human 
retina  and  obtained  a  figure  of  10’  bits/sec.  This  number  seems  high  compared  to  recent 
measurements  of  visual  capacity  [16]  which  indicate  a  minimum  of  50,000  bits/sec  as  the 
required  amount  of  visual  information  necessary  to  discriminate  pictures  for  undersea  work. 

The  measurement  of  capacity  as  bits/sec  should  not  be  confined  with  the  classical  work 
by  Miller  [8]  which  is  concerned  with  capacity  using  units  of  bits/judgment.  Miller’s  law 
showed  that  the  number  of  categories  for  which  simple  objects  of  a  one  dimensional  stimulus 
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could  be  assigned  to  is  consistently  near  7  ±  2  for  a  wide  range  of  different  tasks.  This  means 
approximately  3  bits/judgment  is  the  capacity  of  a  human  for  a  one-dimensional  stimulus. 
Sheridan  [11]  clearly  explains  that  the  human  sensors  act  in  a  multitude  of  dimensions  (each 
dimension  of  which  can  only  distinguish  7  shades  of  stimulus).  Therefore  the  high  capacity 
rates  associated  with  the  reading  of  words  and  the  visual  information  channel  are  due  to  the 
fact  that  these  sensors  are  composed  of  a  multidimensional  channel,  each  dimension  having 
capacity  limited  by  3  bits/judgment. 

These  works  also  support  the  concept  that  when  working  at  capacity,  one  can  increase 
speed  only  at  the  expense  of  accuracy,  and  conversely.  For  example  in  the  study  of  motor 
acquisition  tasks,  Fitts'  law  [17]  concurs  with  the  capacity  measure  of  2-3  bits/judgment  (7 
alternatives)  and  recent  work  [18]  supports  this  fact  in  the  modern  laboratory  environment. 
It  is,  therefore,  desired  in  this  study  to  use  a  variable  of  the  form  I  in  Eq.  (5)  related  to  the 
behavior  model  discussed  previously  to  exhibit  invariance  properties  that  occur  in  informa¬ 
tion  theoretic  models,  as  has  been  discussed  here. 

6.  APPLICATIONS  TO  A  BEHAVIOR  MODEL 

In  order  to  apply  the  previous  concepts  to  the  behavior  model  developed  in  the  last  section, 
it  is  necessary  to  define  a  useful  metric  related  to  tracking  performance  behavior.  If  the 
variable  representing  the  performance  is  a  good  measure,  one  would  expect  ratios  of 
(variations  in  intensity /mean  intensity)  to  show  some  consistency  over  various  levels  of  stress. 

To  achieve  this  goal,  it  is  first  necessary  to  define  the  measure  of  variability  (mean  /i  and 
standard  deviation  <r).  Figure  3  illustrates  the  diagram  that  will  be  used  in  this  procedure. 
A  plot  is  made  of  the  time  derivative  of  e(t)  versus  e(t).  The  densities ffjx),  i  =  1, 2 ...  8,  in 
Fig.  3  are  determined  through  averages  across  replications  of  the  experiment  and  their 
structures  are  not  specified  at  this  point  in  time.  It  is  necessary  to  first  define  the  boundaries 
of  the  e(t)  and  e(t)  axis. 

Definition.  The  phase  plane  boundary  =  fit,  where 


FXA)  =  0.5,  i  =  1, 2 ...  8  (6) 

and  F,(x)  satisfies 

0) 

i  =  1,2 ...  8  and  f/(s)  is  the  appropriate  density  function.  Thus,  the  boundaries  of  the  phase 
plane  are  defined  with  respect  to  median  values  which  can  be  related  to  those  obtained  from 
the  empirical  data  (even  if  the  underlying  densities  are  unknown).  To  determine  the  ±  1  S.D. 
envelopes  about  the  means,  data  are  averaged  across  replications  in  the  following  manner: 

(1)  Find  the  maximum  window  (across  all  replications  k  =  1, 2 . .  ,j). 

(2)  Divide  all  windows  up  into  20  equal  parts  using  0,  x„  x2 . . .  xw. 

(3)  At  each  x„  find  the  F/x,)  value  across  the  j  replications  that  corresponds  to  this  value 
of  x(. 

(4)  Then  average  the  F,  as  follows: 

Fj  =  Z  (Fk)  -  mean  (8) 

*•1 

ft-1 


(9) 
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The  next  step  in  this  procedure  is  to  develop  a  parametric  ratio  which  can  be  used  to  study 
properties  of  the  human  as  an  information  processor.  This  index  should  be  tested  for 
consistency  and  satisfy  properties  similar  to  the  other  variables  discussed.  The  index  chosen 
here  is  specified  by 


I  =  6rjt,  (10) 

where  af.  satisfies  Eq.  (9)  and  fi  is  specified  by  (6). 

7.  EMPIRICAL  VALIDATION  OF  THIS  INDEX 

Two  tests  will  be  conducted  on  the  variable  /  specified  in  Eq.  (10).  First,  a  test  on 
consistency  across  experimental  conditions  will  be  performed  (analogous  to  Weber’s  law). 
Second,  the  calculation  of  processing  bandwidth  will  be  conducted  for  two  experimental 
conditions. 

The  data  base  used  in  this  paper  covers  an  interesting  experiment  in  the  area  of  acceleration 
stress  [12].  As  described  earlier  in  Sec.  2,  human  subjects  were  exposed  to  acceleration  forces 
while  simultaneously  required  to  perform  a  tracking  task.  Figure  5  illustrates  the  stress  levels 
and  durations  of  exposure.  Five  different  tasks  were  chosen  in  this  study  based  on  different 
levels  of  velocity  and  acceleration  profiles  of  the  variable  /(().  A  more  complete  description 
of  the  experiment  can  be  found  in  [12],  For  simplicity  only  the  analysis  details  related  to 
information  theory  are  considered  here. 

The  first  test  of  interest  in  these  data  is  the  level  of  consistency  of  the  variable  I  in  Eq. 
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Table  I.  Values  of  /  (7  subjects  4  replications). 


ff#l 

ff#  2 

ff#  3 

ff#4 

ff#  5 

Exp.  Cond. 

/* 

ff 

/ 

fi 

a 

/ 

fi 

ff  1 

ff  / 

fi 

ff  / 

Static 

0.45 

0.09 

0.2 

0.70 

0.1 

0.143 

1.2 

0.10  0.083 

1.5 

0.11  0.073 

1.6 

0.13  0.081 

Stress 

0.55 

O.U 

0.2 

0.75 

0.12 

0.16 

1.3 

0.13  0.1 

1.7 

0.14  0.082 

1.8 

0.15  0.083 

(10).  If,  indeed,  this  variable  is  a  viable  measure  of  human  tracking,  one  would  expect  a 
Weber's  law  effect  to  hold.  Using  data  from  7  subjects,  4  replications  of  each  experimental 
condition.  Table  1  is  constructed. 

From  Table  1  it  is  seen  that  as  the  task  number  increases,  the  variability  (tr)  increases.  Also, 
going  down  a  column,  both  fi  and  a  increase.  This  is  typical  of  data  from  a  stress  experiment 
of  this  type.  The  ratio  a/fi,  however,  does  not  show  consistency  going  across  a  row  in  either 
the  stress  or  static  condition. 

A  second  calculation  is  conducted  which  involves  processing  bandwidth  using  the  values 
of  Cp  specified  in  Eq.  (1).  If  the  capacity  of  the  human  is  constant,  one  could  write  the 
processing  bandwidth  in  the  following  form: 


^-CyPoftoCl+S/W)],  (11) 

where  i  stands  for  stress  or  static.  The  variable  (1  +  SIN),  related  to  units  of  [1  +  (l//)2],  is 
calculated  and  illustrated  in  Table  2  with  Wt  in  units  of  C,  and  let 

a^ll+d//)2].  02) 

The  data  displayed  in  Tables  1  and  2  are  illustrated  in  Figs.  6a-b.  Table  2  illustrates  some 
important  results.  For  tasks  1  and  2  (ff  #  1  and  ff  #  2),  there  appears  to  be  high  levels  of  Wt. 
For  tasks  3,  4,  and  S  a  different  set  of  the  Wt  values  appear.  This  seems  to  indicate  that  this 
experiment  produces  two  types  of  tracking  behavior  (both  occurring  in  region  I  of  tracking). 
The  high  accuracy  tracking  occurred  for  tasks  1  and  2.  For  the  remaining  3  tasks,  a  different 
type  of  tracking  behavior  is  noted. 

One  hypothesis  to  explain  these  results  found  in  the  data  comes  from  studies  in  muscular 
activity.  Physiologists  speak  of  slow  twitch  muscles  and  fast  twitch  muscles  when  describing 
arm  movements  in  tracking.  The  data  seem  to  reflect  this  type  of  effect.  The  slow  twitch 
muscles  are  used  for  the  accurate  tracking  in  tasks  l  and  2.  The  fast  twitch  muscles  are  used 
in  tracking  where  large  amounts  of  force  are  required.  Since  this  experiment  used  control 
sticks  that  had  large  force  gradients  for  the  higher  numbered  tasks,  this  model  indicates  a 
transition  between  the  two  types  of  muscle  activity  between  tasks  2  and  3.  Thus,  the  modelling 
gives  us  a  clue  as  to  how  to  design  these  tasks  to  investigate  which  type  of  muscle  activity 


Table  2.  W,  m  terms  of  C„  units  and  a. 
ff#l  ff#2  ff#3  ff#4  ff»5 

Exp.  Cond.  a  W,  a  W,  a  V,  a  >  ,  a  W, 
Static  26  0.707  49.9  0.392  145  V.  >87  0.440  133  0.438 


Stress 


26  0.707  40.06  0.624  101  0.498  148  0.461  145  0.463 
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is  affected  by  the  acceleration  stress.  Therefore  by  noting  the  task  numbers  producing 
different  types  of  tracking  responses,  we  can  design  tasks  necessary  to  study  this  effect. 

8.  SUMMARY  AND  CONCLUSIONS 

This  study  discusses  a  behavioral  model  which  incorporates  physiological  stress  and  task 
variables.  Some  similarities  are  noted  between  data  from  this  model  and  information  theory. 
A  proposed  index  for  measuring  human  information  processing  properties  in  the  behavioral 
model  was  tested  for  consistency  and  the  processing  bandwidth  was  calculated  for  two 
experimental  conditions.  Phase  plane  analysis  provided  insight  into  different  types  of  tracking 
which  can  be  used  to  further  investigate  how  humans  use  muscle  activity  in  hand  control. 
By  using  the  behaviorial  model  and  analyzing  the  data,  the  experimenter  developed  a  modus 
operandi  as  to  how  to  design  the  tracking  tasks  for  the  next  experiment  to  study  different 
aspects  of  human  tracking.  Thus  the  use  of  a  mathematical  model  increased  our  knowledge 
about  the  underlying  physical  process  occurring  in  the  data. 

APPENDIX  A 

In  this  appendix  the  closed  loop  error  signal  denoted  as  e{t)  in  the  text  is  analyzed  here 
as  a  stochastic  process  e(t,  x),  where  x  is  the  distribution  of  the  window  variable  and  /  is  the 
time  evolution  of  the  error  signal.  Thus,  with  reference  to  Fig.  3,  if  the  closed  loop  error  signal 
(no  stress)  crosses  the  positive  e(t)  axis  at  time  instant  /, ,  time  can  be  frozen  at  f,  units  and 
the  random  variable  /,(x)  becomes  the  window  distribution  at  that  point  in  time,  i.e., 

/,(x)-e(/„x).  (A.  I) 

Likewise,  if  the  closed  loop  error  signal  crosses  the  positive  e(t)  axis  at  time  t2,  the  random 
variable  /3(x)  is  likewise  specified  by 

f3(x)  =  (d/dt)e(t,x)  |,_,2.  (A. 2) 

To  illustrate  the  difference  in  the  vertical  and  horizontal  axis  in  Fig.  3,  the  variable  z,(t,  x) 
will  denote  the  vertical  axis  and  z2(t,  x)  will  denote  the  horizontal  axis.  Thus, 

z,(t,x)  =  (d/dt)e(t,x)  (A.3a) 

z£t,x)  =  e(t,x)  (A.3b) 

In  this  appendix,  the  assumptions  governing  the  phase  plane  boundaries  will  be  discussed  in 
detail.  With  reference  to  Fig.  3,  the  following  assumptions  will  be  considered. 

Symmetry  Assumptions 

(B.l)  /,(x)  and  /2(x)  are  of  similar  structure  and  symmetrical  about  their  medians,  i.e., 

/.Oh)-/ If-*)  (A.4) 

In  a  similar  manner,  this  property  is  assumed  to  apply  to  /5(x)  and  /6(x). 

(B.2)  /}(x)  and  f,(x)  are  of  similar  structure  and  symmetrical  about  their  medians,  i.e., 

(A.5) 

Similarly,  this  property  is  assumed  to  apply  to/^x)  and  /,(x). 
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Assumptions  governing  the  ellipse 

(B.3)  /3(y)  is  unknown,  but  ffy)  may  not  be  of  the  form  /,(x),  i.e., 

/sO ')#/.(*)  (A-6) 

where  y  =  cx  +  d  is  a  linear  relationship  and  c  and  d  are  constants. 

(B.4)  The  time  evolution  of  the  elliptical  path  is  governed  by 

{z,f/a2  +  (z2flb2=l.  (A.  7) 

(B.5)  In  the  phase  plane,  the  following  dynamical  equation  holds: 

(< d/dt)z2(t ,  x)  =  z,(/,  x)  (A.8) 

(B.6a)  fix),  1  =  1, 2  are  normal, 

(B.6b)  /<x),  1  =  3, 4  are  normal. 

This  appendix  will  examine  the  feasibility  of  assumptions  (B.1)-(B.6)  and  derive  properties 
arising  from  these  conditions.  It  will  first  be  necessary  to  show  the  following  interrelationships 
between  these  assumptions. 

Lemma  1.  (B.6a)  and  (B.4)  together  imply  (B.6b)  is  false. 

Proof.  Solving  (A.  7)  for  z2  yields 


2i~  y/b1  —  (b/a  j*zf .  (A  .9) 

If  z2  has  density  /,(x)  associated  with  this  axis,  then  the  density  of  z,  subject  to  (A.9)  can  be 
determined  by  coordinate  transformation  [13]: 

f(y)  =/,[(* 2  -  (blafzmd/dzMb2  -  (b/afz^'%,  (A.  10) 

or 

f(y)  =  (b 2l(^/lna 2a)(y/(b2  -  (b/afy2fm) (A.l  1) 

If  one  now  attempts  to  write  this  in  the  form 

f3(y)  =  KMj2x<ry)e~<r-'>W),  (A.12) 

where  y  is  a  linear  function  of  the  variable  x,  one  sees  this  cannot  be  done  due  to  the  term 
y/(b2  -  ( b/afy)om  in  (A.  1 1).  Thus  the  density  f(y)  derived  by  assumptions  (B.6a)  and  (B.4) 
imply  that  (B.6b)  is  not  satisfied.  Next,  assumptions  (B.4)  and  (B.S)  will  be  examined  to  see 
if  they  give  rise  to  some  level  of  consistency. 

Lemma  2.  The  assumptions  (B.4)  and  (B.S)  are  consistent  with  /((x)  arbitrary  and  /3(x) 
satisfying  (A.  10). 

Proof.  Combining  Eqs.  (A.7)  and  (A.8)  yields 


(d/dt)z^t,  x)  =  s/a2-(albfz% 


(A.  13) 
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(A.  14) 

(A.  15) 

(A.  16) 


ZiO,  Mi)  =  a  cos(a(t  -  t0)/b) .  (A.  1 7) 

Thus  the  stochastic  process  may  evolve  in  time  via  the  parametric  elliptical  equations 
(A.  16- 17).  It  is  noted  at  this  point  that  the  density  functions /,(x)  and  f}(y)  are  still  arbitrary. 
From  the  derivation  in  Lemma  1,  however,  it  is  noted  that  if_/i(x)  is  specified,  then /3(y)  has 
a  dependence  through  the  coordinate  transformation  [Eq.  (A.  10)]. 

In  summary,  the  following  assumptions  concerning  the  stochastic  process  e(t,  x)  can  be 
stated  which  are  consistent  with  all  the  analysis  conducted  in  this  paper: 

(I)  In  the  phase  plane  axis,  the  following  relationship  holds  for  all  x: 

(dldt)(e(t,  x))  =  (d/dt)z^t,  x)  =  z,(f,  x) .  (A.  18) 

(II)  The  time  evolution  of  the  elliptical  path  is  governed  by 

(z\)Ha*)  +  (zl)/(b*)=\.  (A.  19) 

(III)  The  density  of /,(x)  may  be  specified.  From  [12]  it  is  known  that  this  function  should 
not  be  assumed  to  be  normal.  If /,(x)  is  specified  then  the  density  f,(y)  must  satisfy 

Uy)  -  (b/afzftymz,  )•[((**  -  {b/aYzm  .,1,.^  (A.20) 
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